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ON SUM CRITERION FOR THE EXISTENCE OF 
SOME ITERATIVE METHODS 


K. RAUF, O. T. WAHAB, J. O. OMOLEHIN, I. ABDULLAHI, O. R. ZUBAIR, 
S. M. ALATA, I. F. USAMOT AND A. O. SANUSI 


ABSTRACT. In this paper, we evoke an existence for two iterative 
methods of the system of linear equations. Our results agree with 
existing results and suggest a strong procedural condition for the 
convergence of the system of linear equations. 


1. INTRODUCTION 


We consider the system of linear equation 


(1) Ar =b 
where 
Qi, Q12 ... Qin Tı bi 
Q21 Q2 ... Gan T2 be 
A= a ; , b= 
Ani an2 ++» Ann Tn bn 


Problems concerning system of linear equations are often occurred in 
engineering and sciences. Many researchers have developed several 
methods for solving the problems. Noor et al. [2] had used decom- 
position method for solving the system (1). Babolian et al.[6] have 
used Adomian decomposition method to derive an iterative method 
which is similar to the Jacobi iteration. Allahviranloo [4] used Ado- 
mian decomposition method for fuzzy system. Keramati [7] and Liu 
[8] have used homotopy perturbation method to derive some iterative 
methods for solving the system. Improving Newton-Raphson method 
for non linear equations by modified Adomian decomposition method 
was discussed in [3]. In this paper, we evoke a strong procedural con- 
dition for the decomposition method in [2] which is sufficient for the 
existence of system (1). 
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2. PRELIMINARIES 


We express system (1) as: 


Ly =l E a1 )£1 — QAi2ť2 — ... — AlnTn + bi 

LQ = — Q21£1 + (1 — az22)£2 — . . . — AmTn + b2 
(2) 

Ln = — Ani®1 — AneT2 — ... + (1 — ann)En + bn 


If we set fij = Oi; — aij; 1 < i,j < n, where 46;; is the kronecker delta 
|]. Equation (2) can be written as 


£1 = Pati Pete — Pia + bı 
Lo = —P21£1 + Boo — . . . — Bon Ln + b2 
(3) i 
Ln = — bni — Br2X2 Se AP Pinte + bn 
This implies 
j=1 


And in matrix form as 
(5) x= Bxr+b 
where B = [By], 1<i,j <n 


Definition 2.1. Let (R",a(-)) be a metric space. Let T be a 


mapping such that T : R” — R”. T is called a contraction on R” if 
there is a positive real number a < 1 such that for all x,y € R” 


d(Tx,Ty) < ad(z,y). 


The point x € R” is fixed if it coincides with its image. i.e x = Tx. 
Now, define a mapping T on (5) which map R” into itself, then 


(6) Tx = Bz +b 


and we can write 


(7) Tx; = )_ Biytj + bj, 1=1,2,...,n. 
j=l 
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2.1. Decomposed System. Let h 4 0 be a parameter, for any s- 
plitting matrix Q and an auxiliary matrix H, we can decompose the 
system (1) as 


(8) Oi GHAAOPEED 
Let Q = D = diag(aj;), i = 1,2,..., then we can compare (8) with (5) 
if we let 
B' = -D (hH A — D) and b = D"'hHb 
so that 
(9) aD — Biol 4y 


where m is a positive integer representing the number of iterations. 
This is called the Adomian Jacobi method [5]. 

If we let Q = -L + D 

where L is the lower triangular matrix with principal diagonal all zero. 
From (8), we have 


af) — (I —h(D — L) HA)z™ + h(D — L) Hb 


(10) gmt) = B'r i p" 


where B” = (I — h(D — L)~'HA) and b” = h(D — L)~'Hb 
This is called the Adomian Gauss-Seidel method [5]. 


3. MAIN RESULTS 


In this section, we obtain the generalization of sum criterion and 
deduce some results which are sufficient for the existence of (1). 
The following Lemma is significant to our main results. 


Lemma 3.1. Let R” be a metric space with metric dp on R” and T 
be a mapping of R” into itself. Then, T is a contraction if 


n 


(11) 5 S Bal? <1, 1<p<o 


i=1 j=l 


Proof 
Let R” be a metric space, for x,y € R”, the metric d, on R” is defined 
by 


1 


(12) d,(x,y) = > |£; -u) Terso 
i=1 
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Let T : R” — R” and x,y € R”. Using (6), we have 


1 


(>: Ta- TuP) 
> Bilas Yj) ) 
< E La De yj) 


i=1 j=l 


d (Tz, Ty) 





e iM 





aa 











er 











A 1 
P 


= bss ae ar hs [£j — wr) (Schwarz’s inequality) 
i=1 j=1 p 
o [isl T dp(x, y) 
i=1 j=1 


Choose a = 2s Xa Bal) j 
hence, 

d,(T x, Ty) < ad,(x, y) 
Therefore, the mapping T on R” is a contraction. 


Remark. If the condition (11) is applied on either (9) or (10), then we 
obtain 


(13) > leal? < lox? 

i=1 j=1 
Inequality (13) implies that the elements in the principal diagonal of 
matrix A must be large, and it will be called p-th sum criterion. 
The condition (11) is sufficient for the convergence of (6) while (13) is 
sufficient for the convergence of (9) and (10) but the manner at which 
they converge differ depending on the choice of p. 
We consider the following cases 


Case I: Suppose p = 1 in equation (13), we obtain 


n 
` laij| < Jal 
i=1 
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for 7 = 1,2,...,n and it is called column sum criterion 
Case II: Suppose p = 2 in equation (13), we have the euclidean metric 
and the condition is called the square sum criterion given by 


n n 
2 2 
5 ` Qij < Qü 


i=1 j=1 
See [1] for similar work. 


Case II: When p = 3 with the metric defined on R”. For (13), we obtain 
the condition 


n n 
3 3 
` ` la|” < jaul 
i=1 j=1 
which is called the cubic sum criterion 


Case IV: If p = œo with the metric dœ on R”, we obtain the condition 
sup Š. laij| < Jaxl 
a j=l 
This is called the row sum criterion. 


We justify the aforementioned cases with the following examples: 
Example 1: 
Consider the system of linear equation 
4a —12- e = 0 
—X1 + 4x5 — T3 — T5 =5 


—T2 T 473 — = 0 





=b T 404 e 6 








—Tə — T4 + 4T5 — T6 = —2 
—g£3 — £5 + 4x6 = 6 
Solution 

6 

j=1, dle = laa + |azi| + |aa1| + |asi| 4 laal = 7 
3 
lail > 7 for i = 1,2, 3,4,5,6. 

Case I is satisfied. 
6 6 

2, 2, 2, 2, 2,2, 2,2,2, 2,2,2,2 _ 1 





T 
|axi| > 9? for i = 12: BA 5.6. 
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Case II is satisfied 


Boa 


ease: 
Qij A127 a1 T023 T425 A32 TA36 T41 TA45 A521 A54 A56 T63 G64 = 39 





i=1 j=1 


7 
lail > 39° for i = 1,2,3, 4,5,6. 
Case III is satisfied 


6 
1 
sup X` laij] = n 1 
(3 j=1 


1 , 
|as| > T for i = 1,2,3, 4,5,6. 
Case IV is satisfied 


Example 2. 


Consider the system 
2x1 + Lo T T3 =4 








£i 2X9 v3 = 4 





U1 +%2+24%3=4 


Solution 
Using the cases above, we have 
3 
S- |aij] = 1 = Jail, for ¢= 1,2, 3. 


j=l 


This gives a non-expansive type. 


5 s3 6 
DIED E E E E A O ES Pn 
5 5 Qij = Qiz T A713 T Q21 T Q23 T Q31 T a32 = 7 


i=1 j=l 





6 
lail < T for i = 1,2,3. 


Case II is not satisfied. 


3 3 3 
Bh ma Be ge E: E he EE 
` 5 Qij = Qiz T 413 T Q21 T Q23 T Q31 T wg 


i=1 j=l 





3 
lai] > T for i = 1,2,3. 


i 1 
sup Š |aij| = 5 
{A j=1 


Case III is satisfied. 
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1 
lail > z for i = 1,2,3. 


Case IV is also satisfied. 
Example 3: 
Consider the system 

Tı + T3 = 2 


—-% T T2 = 0 





£i T 2X9 — 323 = 0 


Solution : 


4 
j=1, $ layl = lan] + lasi] = 3 


i=l 
4 
|axi| < 3? for i = 1,2,3. 


Case I is not satisfied. 





23 
lal < D for i = 1,2,3. 


Case II is not satisfied 


5-3 
47 
Bo 2 nies is ES ETEN oe i a a 
5 X dij = ala + afs + 03) + a33 + a3, + O59 = zy 


i=1 j=l 





47 
laul < aR for i = 1,2, 3. 


Case III is also not satisfied 
3 


sup ò | laij| = 1 
(3 j=1 
lail = 1, for i = 1, 2,3. 


This is non-expansive. 


3.1. Conclusion. It is generally known that if the system (5) of n 
linear equations satisfies at least one of the cases I, I, II and IV, then, 
it has precisely one solution. 

It would be observed that the example 1 is satisfied for all the cases, 
and hence, both the Jacobi and Gauss-seidel iterations converge. In 
example 2, case I gives a non-expansive while case II is not satisfied, 
and therefore, the Jacobi iteration is not converge but Gauss-seidel 
iteration does. Only case IV gives a non-expansive in example 3, while 
cases I, I] and III were not satisfied, so, Jacobi iteration converges while 
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Gauss-seidel iteration diverges. On this note, the Jacobi and the Gauss- 
seidel iteration are not comparable, though the Gauss-seidel iteration 
converges faster than the Jacobi iteration but we are concerned on the 
criterion at which the iterative methods converge. 
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